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Abstract: 



In this note we solve a proljleni about the rational representability of hypergeometric terms which 
represent hypergeometric sums. This problem was proposed by Koornwinder in [4]. 

1 Hypergeometric Functions and Zeilberger's Algorithm 

Zeilberger's algorithm ([6]-[7], see also [3], [1]) determines recurrence equations for hypergeo- 
metric functions 



whose upper parameters and lower parameters Pi~ are rational-linear in a variable n, when- 
ever the term ratio 



is a rational function in both n and k. As usual {a)k = a{a + 1) ■ ■ ■ {a + k — I) denotes the 
Pochhammer symbol. We call the summand A^x^ a hypergeometric term. The resulting 
recurrence equation has polynomial coefficients with respect to n. If it is of first order, the 
sum has a rational term ratio with respect to n, and hence is itself a hypergeometric term. 
In [4] Koornwider asked the question whether an application of Zeilberger's algorithm might 
generate a hypergeometric term whose upper and lower parameters are not rational assuming 
the parameters of the input summand are rational: 

Problem 6.1. If Zeilberger's algorithm succeeds, can S{n)/ S{n—1) then always be 
factorized as a quotient of products of linear forms over Z in n and the parameters? 

In this note, we will answer Koornwinder's question in the negative, by providing several 
counterexamples. 

Note that Koornwinder's question in principle is independent of Zeilberger's algorithm, and 
asks whether there are hypergeometric sums that can be represented by hypergeometric terms 
with nonrational parameters. None example of this type can be found in the literature, see in 
particular the rather extensive mathematical dictionary on hypergeometric function identities 
[5]. Nevertheless we will use Zeilberger's algorithm to find our counterexamples. 
In [1], we presented a Maple implementation for an extension of Zeilberger's algorithm, which 
is available through the sumtools package of Maple V.4 (or through the share library package 
summation). In this note we use a new implementation, developed in the book [2], which can be 
obtained from http://www.zib.de/koepf/code and http://www.zib.de/koepf/simpcomb. 
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but the recurrence equation calculations can also be done with the sumtools package of Maple 
V.4. 

In [1], we gave the following example: 
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This example shows that despite the radicals of the input series, the output is rational. 
More examples of this type can be found in [5]. The above example can be generated 
by the Maple command (with the package code, or in Maple V.4, after with(sumtools) : 
readlib ( ' sum/ simpcomb ' ) : ) 

> sumrecursiondiypertermC [1/2, -2*n] , [2*n+3/2] ,3+2*sqrt(2) ,k) ,k,S(n)) ; 

-(8n + 9)(8n + 7)S(n + l)+4S(n)(4n + 5)(4n + 3)=0 

or by 

> sumrecursionC (-1) "k*binomial(2*n,k) *binomial(2*n+k+l ,k)/ 

> binomial (4*n+2*k+2,2*k)*(3+2*sqrt (2) )~k,k,S(n)) ; 

-(8n + 9)(8n + 7)S(n + l)+4S(n)(4n + 5)(4n + 3) = 0. 

Prom the resulting recurrence equation which is valid for the above sum, one can easily read 
off the term ratio S(n + 1) /S(n), hence the parameters of the hypergeometric sum. Using the 
package code, this can be automatically done by 

> Closedf orm(hyperterm( [l/2,-2*n] , [2*ii+3/2] ,3+2*sqrt(2) ,k) ,k,n) ; 
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But Koornwinder's question is a different one: He asks whether rational input (the summand 
of a hypergeometric sum) can generate nonrational output (a hypergeometric term represent- 
ing the sum under consideration). Obviously rational input generates a rational recurrence 
equation. But the question remains whether the term ratio S{n + 1)/S{n) equivalent to this 
recurrence equation can be factorized rationally. All known examples in the literature are of 
this type. Note that rational factorization can be done algorithmically and is accessible in 
computer algebra systems like Maple. 

2 Hypergeometric Product Formulas 

We have discovered our counterexamples in connection with the examination of product for- 
mulas. 

As an example we consider Clausen 's formula 



2^1 



a, b 
a + b+1/2 



3^2 



2a, 26, a + b 
a + 6 + 1/2, 2a + 26 



(1) 



Normally Clausen's formula is deduced by showing that both sides satisfy the same differential 
equation (of third order) with respect to x. Having Zeilberger's algorithm at hand, there is 
a second possibility, though. Represent the left hand side as Cauchy product. This is a 
representation by a double sum. Find a recurrence equation for the inner sum with respect 
to the summation variable of the outer sum. If the resulting recurrence equation is of first 



2 



order, then it determines a hypergeometric term, the summand of the right hand sum. Hence 
the right hand sum constitutes a hypergeometric function whose parameters can be read off 
directly from the detected recurrence equation. 

For the coefficient of of the Clausen product on the left hand side of (1), represented as 
Cauchy product, we get 

> sumrecursionC 

> hypertermC [a,b] , [a+b+1/2] , 1 , j ) *hyperterin( [a,b] , [a+b+1/2] ,l,k-j) , j ,S(k)) ; 

-(A; + l)(2a + A; + 26)(2a + 26+l + 2A;)S(A; + l) 
+ 2S(A;)(26 + A;)(A; + 2a)(a + A; + 6)=0 

This recurrence equation is valid for the coefficient of the outer sum, and hence this coefficient 
is given by 

> Closedf ormC 

> hypertermC [a, b] , [a+b+1/2] , 1 , j ) *hyperterin( [a,b] , [a+b+1/2] ,l,k-j) , j ,k) ; 



Hyperterm ^[2 6, 2a,a + 6] 



a + 6 + -,2a + 26 



This hypergeometric term is the coefficient of x^ of the right hand sum of (1), which finishes 
the proof of (1). 

All the usual product theorems of hypergeometric series can be deduced by this method, see 
[2], Chapter 7. Note that the given procedure has the advantage over the differential equation 
approach that it generates the right hand sides, given the left hand products. We will consider 
more examples soon where we detect such representations without prior knowledge. 
The calculations 

> Closedf orm (hypertermC [] , [] jX, j ) *hyperterm( [] , [] ,y ,k-j ) , j ,k) ; 

Hyperterm([],[],a; + ?/,A;) 



and 



> Closedf orm(hyperterm( [a] , [b] ,x, j)*hyperterm( [] , [] ,-x,k-j) , j ,k) ; 

Hyperterm( [ —a + 6], [ 6], — k ) 



generate the product identities 



and 
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X 
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(2) 
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Obviously the first one is the addition theorem of the exponential function, and the second 
one is Kummer's identity. 

The computations above deduced Equations (2) and (3) by converting 



k-j 



j=0 



j! {k-jV. 



to {x + y)^ /k\, and by converting 



^ {a), x^ {-x)^-i 
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to 



{b-a)k i-x)^ 



{b)k k\ ' 

respectively. 
The calculation 

> sumrecursion(hyperterm( [a] , [b] ,x, j ) +hyperterm( [a] , [b] ,-x,k-j ) , j ,k) ; 
{k + 2) {k + 1 + b) {k + b) {k + 2b)S{k + 2) + S{k) {2a + k) {2a - k - 2b) = 



shows that 
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Since the product is an even function, this can be also computed by 

> Closedf orin(hyperterm( [a] , [b] ,x, j)*hyperterm( [a] , [b] ,-x,2*k-j) , j ,k) ; 

1 



Hyperterm ( [ —a + b,a 



3 Counterexamples 

Now we investigate examples with nonrational results. A simple class is given by 
> Closedf orm (hyperterm ( [a] , [a-1] , 1 , j ) *hyperterin( [a] , [a-1] , 1 ,k-j ) , j ,k) ; 



Hyperterm y 



2a-^-^V9-Sa,2a-^ + ^V9-Sa 



1 



2a- ---V9-8a,2a- - + - V9-8a 



,2,k 



in particular for a = 3/4 

> ClosedformC 

> subs (a=3/4, hyperterm ( [a] , [a-1] , 1 , j ) *hyperterm( [a] , [a-1] , 1 ,k-j ) ) , j ,k) ; 



Hyperterm I 



VL 2 



,2, A; 



generating the identities 



1^1 



2-f 2 



2a - 3/2 



,2a- 3/2 + 



2a - 5/2 - 2a - 5/2 + 



2x 



and 



1^1 



3/4 
-1/4 



2 ' 2 
_ 1 _V3 _ I 



2x 



2^2 1 ^ 

2 ^' 2 

respectively. The latter gives a simple counterexample for Koornwinder's question for which 
radicals cannot be avoided. In terms of summations, it reads as 



E 



(3/4), 



(3/4),- 



i-ik'{3/i)j{5/4-k)j{-k)j {-ly 



j=0 

l-Ak 



(-1/4),) j! (-l/4),_, (^--j)! 

' 3/4, 5/4 -A;, -A,- 



A;! 



E 



A;! 



3-P2 



1/4- A;, -1/4 



-1 



(1/4 -A;),- (-1/4),- 

2k 



(V3/2), (-a/3/2). 



(a/3/2 - 1), (-V3/2 - 1), A;! ' 
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and the generic example involving the parameter a reads as 

(«)j Wfc-j ^ a+k-l ^ {a)j {2-a-k)j {-k)j j-iy ^ 

^^{a-l)j)j\ {a-l)k-j{k-j)\ {a-l)k\f^^ - l)^ {l-a-k)^ j\ ^' 

"7 " (-5/2+2a-^%^J^ |_5/2+2a+^^j^ k\ ' 




(a-1) /c! 



These are 3F2-evaluations by nonrational hypergeometric terms. 

Note that the hypergeometric representation given in (4) can be deduced by the Maple com- 
mand 

> simpcombCSumtohyper (hypertermC [a] , [a-1] , 1 , j ) *liyperterm( [a] , [a-1] , 1 ,k-j ) , j ) ) ; 

{a + k — 1) Hypergeom( [a, —a + 2 — k, —k], [a — 1, —a — A; + 1], — 1 ) 

(a-l)r(A; + l) 

Similarly for n = 2, 3, . . . the products 

\ 2 




lead to a recurrence equation of first order with respect to k with polynomial coefficients of 
degree 2n in k that have no proper factorization over Q, e.g. for n = 2 

> sumrecursion (hypertermC [a] , [a-2] , 1 , j ) *hyperterm( [a] , [a-2] , 1 ,k-j ) , j ,S(k) ) ; 

-{k + l)(k'^ - Uk^ + 8ak^ + 67 k"^ - 80ak^ + 24a^ k^ - nSk + 232ak + 32a^ k 
- 152 a'^k + 208 - 96 + 16 - 192 a + 64)S( A; + 1 ) + 2S( A; )(-22 k 
+ 96 a A; - 56 a A;2 - 104 A; + 80 + 31 - 32 a - 10 k^ - 64 + 16 
+ 32a^ k + 24a'^ k'^ + k^ + 8ak^) = 

or for n = 3 

> sumrecursion (hypertermC [a] , [a-3] , 1 , j ) *hyperterm( [a] , [a-3] , 1 ,k-j ) , j ,S(k) ) ; 

-{k + l) (-264 k^a- 8448 a - 5136 k + 289 + 6180 k^ + 2180 ak^ - 9216 

- 2016 a:' k^ - 1032 k^ + 7808 a'^ k + 3712 - 27 k^ + 2304 + 12 A;^ a + A;^ 

- 1968 a'^k- 768 + 13984 ak + 4102 k'^ + 192 k + 240 k'^ + 160 k^ 
+ 60 k^ - 15000 a'^k + 12352 - 8232 ak'^ - 1533 A;^ + 64 a^)S( A; + 1 ) + 2 
S( A; ) (-204 A;^ a - 768 a - 504 k + 169 k^ + 3444 k^ + 1244 ak^ - 3264 

- 1536 k^ - 792 + 4256 A + 1984 - 21 A^ + 12 A;'^ a + A*^ - 1488 k 

- 576 a''^ + 3064 a A; + 982 A;^ + 192 k + 240 k^ + 160 k^ + 60 k^ 

- 5496 A; + 2560 - 3156 a A^ - 627 A^ + 64 a^) = 

To present some other results of similar type, we calculate the product 



-X 
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which is representable by a hypergeometric function with parameters using square roots, but 
in a rather complicated way: 

> sumrecursion( 

> hypertermC [a, a+b] , [a-2] ,x, j ) *hyperterm( [a,a+b] , [a-2] , -x,2*k-j ) , j , S(k) ) ; 

-{k + l){-7ak'^ + 5k'^b + k'^b'^ + a'^k'^ + 8k'^ -2abk'^ + 21ak + 4:abk-12k 

-2ka'^b-13a'^k + 2ka^ -5kb-kb'^ + 4:-12a + a^ -Qa^ + 13a'^)S{k + l) 
+ S{k){b + k + a){a^ -4:a^ + 2ka^ -lla'^k + a^k'^ -2a'^b + a'^ -2ka'^b 
+ 2a-lak'^ + lak + 2ab-2abk^ + %k^ + 'akb + k^b^ + kb^ + hk'^b + 4:k) 
= 

generating the hypergeometric term 

> factor (ClosedformC 

> hypertermC [a, a+b] , [a-2] ,x, j ) *hyperterm( [a, a+b] , [a-2] , -x,2*k-j ) , j ,k) ) ; 



Hyperterm 



, I2a^ -lla"^ -2a'^b + 7a + hb + l)^ + 4.-8 
a + b,- 



2 A 
1 2a^ -11 a'^ -2a'^b + 7a + 5b + b'^ + 4: + B 



2 A 

1 -ISa'^ + 4:ab + 21 a - 12- 2a'^b + 2a^ -b'^ -5b- B 

2 A ' 
1 -13a^ + 4:ab + 21a-12-2a'^b + 2a^ -b^ -5b + B 



with 



A 



A ■= b'^ + 8 + 5b - 2ab + a'^ - 7 a 



and 

B := ({a + l + b){a^ + %a^ + 4:a^b^ + 150^ b-2Aa-Mab-%ab^ +24:b + %b^ + b^ + 1&)) 

Similarly 

> sumrecursionC 

> hypertermC [a, a+b] , [a-3] , 1 , j ) *hyperterm( [a, a+b] , [a-3] , -1 ,2*k-j ) , j , S(k) ) ; 

-(A; + l)(36 - 132a- 132A; + 58a^ -84a6A;2 - UAa^ + - 5Ak b + ^abk 
+ 168 k^ + ?,ab^ + ?,a^ k^ b^ + 193 - %k^ b^ + A2a^ bk^ + 24:abk^ 
+ 18a^bk-^a'^bk^ -Qa^bk^ -Iha'^k^ -A2a^ k^ -i%a^k + a^ k^ + 2,a^ k^ 
+ 56 a A;^ + 144 k^ - 351 + Blfc^ 5 + 3 A,- - 3 a'^ 6 - 12 a'^ - 2 b^ 

- 255 ak^ + llQka^ + 2,kH^ -^ak^b^ -ia^ k b^ - k^ b^ - ASk^ + 6ab^ k 

- Ibkb"^ + 24k'^ b"^ - 49kb + 337 ak - 32k^ b)S{k + 1) + S{k) {b + k + a){6a 
+ 12k + 22a'^ -Uabk'^ -15a^ + 21ka'^b-30abk + 123a'^k'^ 

+ 3ab'^ k^ + 3a'^ k^ b^ -ha^ -9 k^ b^ + 33a^ bk^ + 2Aabk^ + &a^ bk 
-3a'^bk^ -Qa^bk^ -15a^k^ -39a^k^ -33a'^k + a^k^ + 3a^k'^ + 5^ak^ 
-15a^b-ma^k-15k^b + 3a^k-3a^bk-9a^ +kb^ + &ab-87ak^ 
+ 89ka^ + 3a'^ kb"^ - k^b^ - 48k^ - 3a^b-3ab'^ k + 6kb^ - 3k'^b'^ + 17kb 
-5ak-32k^b+12a^ b) = , 



1/2 
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and first order recurrence equations, more and more complicated with increasing n, can be 
determined iteratively for the coefficients of 



2^1 



a,a + b 



a — n 



a,a + b 



a — n 



We finish this note with a rather simple family of examples. If n G N is a nonnegative integer, 
then the polynomials {k G N) 

/ a + n, —k 



2^1 



satisfy a recurrence equation of first order with respect to k with polynomial coefficients of 
degree n in A; that have no proper factorization over Q. As a particular case one has for n = 3 



> suinrecursion(hyperterin( [a+3 , -k] , [a] ,x, j) , j ,S(k)) ; 

{-2a-6x'^ ka^ + 6x^ ka + 3x^ ka^ + 2x^ k + 2x^ a-3x'^ + + 3x^ k'^ 

— Ibx'^ka — ix'^k'^a + ix^k'^a — 6x'^k'^ + ix^a^ + x^a^ + 9xka + 6xa + 6xk 
-Sa'^ - + 3xka'^ + 3xa^ + 9xa'^ -Qx'^ a-6x'^ k -9x'^ a'^)S{k + 1 ) + S( A; ) 
{x -l ){-2a + 6x + 6x^ -6x^ ka^ + 12x^ ka + 3x^ ka"^ + 11 x^ k + 11 x^ a 
-Sx'^a^ + x^k^ + 6x^k'^ -21x'^ka-3x'^k'^a + 3x^k'^a-6x'^k^ + 6x^0^ 
+ x^a^ + 9xka + 16xa + 6xk-3a'^ -a^ + 3xka'^ -12x'^ + 3xa^ + 12xa^ 
-24x'^a-18x'^k-15x'^ a^) = 
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